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SECTION  I 


{ 


INTRODUCTION 


Despiue  the  considerable  analytic  and  experimental  work  on  compressiole 
turbulent  boundary  layers  vdilch  has  been  performed  In  recent  years,  a  ratls- 
factory  and  general  description  Is  still  lacking.  The  basic  complexify  of 
turbulent  flow  has  led  those  seeking  to  describe  the  compressible  cas^  to 
hypothesize  similarities  between  the  compressible  boundary  layer  and  the 
constant  density  boundary  layer,  where  at  least  a  larger  body  of  experimental 
measurements  Is  available. 

No  attempt  will  be  made  here  to  I'evlew  the  numerous  analytic  lnvestlga> 
tlons  of  compressible  turbulent  flow  which  have  been  made  in  recent  years; 
however,  mention  of  some  recent  references  for  background  In  order.  The 
work  by  Spalding  and  Chl(^)  Included  a  rather  comprehensive  review  of  current 
theoretical  methods  as  well  as  a  seml>empirical  theory.  The  transformation 
approach  developed  by  Coles^^)  has  been  modified  by  Crocco^^)  and  extended  and 
applied  by  Barontl  and  Libby .  EconomosC^)  recently  developed  a  phenomeno¬ 
logical  treatment  for  solving  the  compressible  turbulent  boundary  layer  ylth 
mass  addition  which  Is  a  generalization  of  the  work  of  Barontl  and  Llbby^^' . 

Another  approach  to  this  complex  problem  Is  to  adapt  methods  which  have 
proven  successful  In  solving  the  laminar  boundary  layer.  In  this  approach. 
Investigators  have  used  an  eddy  viscosity  concept  to  express  the  Reynolds 
stress  terms  which  occur  in  equations  governing  turbulent  boundary  layer  flow. 
The  result  is  a  set  of  equations  for  the  turbulent  boundary  layer  which  differ 
from  the  equations  governing  laminar  boundary  layer  flow  solely  in  that  mole¬ 
cular  viscosity  and  Prandtl  number  are  replaced  by  an  "effv^ctlve"  viscosity 
and  Prandtl  number.  If  one  can  determine  empirically  what  values  to  assign 
the  effective  viscosity  and  prandtl  nunber,  it  becomes  possible  to  solve  the 
mean  characteristics  of  the  turbulent  boundary  layer  using  well  established 
methods  for  the  solution  of  the  laminar  boundary  layer.  Eddy  viscosity  models 
for  wall  boundary  layers  have  been  developed  by  Klelnsteln'”! ,  Smith  and 
Cebecl(^),  Patankar  and  Spalding(^),  and  Sontow6kl(^)  to  name  some  of  the 
more  recent  contributors.  Each  of  these  models  was  developed  for  the  zero 
pressure  gradient,  zero  mass  injection  case.  All  were  shown  by  Martelluccl, 
Rie,  and  Sontowski'^^'  to  present  reasonable  agreement  with  data  even  when 
moderate  rates  of  mass  injection  were  considered.  Cebeci,  Smith,  and 
Mosinkis^^^^  recently  modified  the  viscosity  model  of  Ref.  7  to  include  the 
effect  of  pressure  gradient. 

Attempts  to  compare  the  theoretical  work  with  compressible  turbulent 
flow  experiments  have  been  limited  not  only  by  the  paucity  of  detailed  data 
but  also  by  the  incomplete  experimental  definition  of  the  mean  profiles  across 
the  boundary  layer  in  high  speed  compressible  flow.  This  Is  especially  true 
for  the  hypersonic  case  and  where  mass  injection  and  other  effects  such  as 
temperature  junps,  tdilch  cause  non  similar  profiles,  are  present. 


Because  of  the  complexity  of  the  more  rigorous  analytic  formulations,  the 
numerical  difficulties  encountered  in  the  solution  of  ablating  vehicle  boundary 
layers  (due  to  the  ablation  rates,  discontinuities,  etc.),  and  the  excessive 
computer  times  required  to  generate  a  solution,  vehicle  designers  are  forced 
to  use  engineering  analyses  which  are  formulated  to  solve  the  boundary  layer 
equations.  Due  to  the  uncertainties  introduced  by  these  empirical  relations, 
a  high  degree  of  rigor  in  the  mathematical  solution  of  the  equation  is  not 
warranted.  As  a  consequence,  for  the  analysis  of  a  turbulent  boundary  layer 
with  mass  injection,  most  authors^^^”^^)  have  utilized  additional  assumptions 
which  have  greatly  simplified  the  mathematical  solution  of  the  equations. 
However  agreement  between  existing  theory  and  experimental  data  for  a  wide 
range  of  conditions  is  poor^l^"!^!.  A  need  exists  for  an  analysis  that  mini¬ 
mizes  the  dependency  on  the  assumptions  usually  employed  to  simplify  the 
solution  of  the  equations,  but  which  retains  the  essence  of  the  more  rigorous 
formulations.  These  engineering  analyses  were  generated  by  the  various  aero¬ 
space  companies' 20-25)  ^he  prime  purpose  of  computing  the  vehicle  loads 
and  drag  (which  implies  an  accurate  definition  of  the  local  properties). 

Since  these  empirical  analyses  were  formulated,  additional  data  have  become 
available  which  can  be  used  to  validate  and/or  update  some  of  the  empiricism 
that  is  Inherent  in  these  methods. 

•The  objective  of  the  current  effort  is  to  re-evaluate  the  influential 
mechanisms  of  interaction  between  re-entry  body  heat  shield  ablation  and 
aerodynamic  parameters,  based  upon  data  obtained  under  the  preceding  STREET-G 
contract,  other  ground  test  data,  with  analytic  results,  and  with  applicable 
flight  test  data.  The  specific  objectives  of  this  research  are  to  improve 
(1)  the  empirical  correlations  which  account  for  mass  injection  effects  on 
the  viscous  layer  thicknesses,  (2)  the  basic  turbulent  relations  for  pre¬ 
dicting  the  viscous  properties,  and  (3)  the  validity  of  the  Crocco  integral 
for  turbulent  flows  with  mass  injection.  The  influence  of  the  improved  corre¬ 
lations  on  the  vehicle  drag  for  some  typical  flight  vehicles  1*  also 
ascertained. 


SECTION  II 

EXPERIMENT  DESCRIPTION 


Since  the  primary  data  that  are  to  be  used  in  the  current  research  were 
obtained  and  reported  under  the  preceding  STREET-G  program^^S)^  only  a  brief 
description  of  the  test  facility,  the  test  model,  and  the  boundary  layer 
probes  is  in  order. 

2.1  TEST  FACILITY 

The  test  program  was  conducted  in  the  Arnolc’  Engineering  Development 
Center  (AEDC)  Tunnel  B  at  Mach  8.  This  facility  Is  a  continuous  flow,  hyper¬ 
sonic  wind  tunnel  with  a  50-inch  diameter  test  section.  Interchangeable 
axisymmetric  contoured  nozzles  provide  flow  at  either  Mach  6  or  8.  The 
tunnel  is  a  closed  circuit  type  which  operates  over  a  range  of  pressures  with 
air  supplied  by  a  central  compressor  system.  A  modal  injection  system  pro¬ 
vides  the  capability  for  quick  model  changes  without  Interrupting  tunnel 
operation.  The  tunnel  flow  conditions  provide  free  stream  Reynolds  numbers 
from  0.30  x  10^  to  3.8  x  10^  per  foot  at  Mach  8  with  a  maximum  total  tempera¬ 
ture  of  1350^.  A  sketch  of  the  tunnel  layout  and  a  summary  chart  of  the 
performance  characteristics  of  the  facility  are  shown  in  Figure  1. 

Model  support  and  sting  assemblies  provide  a  30-degree  pitch  range 
capability.  A  ntinber  of  conventional  internal  strain  gage  balances  is 
available  at  the  facility  to  measure  data  in  mai^  ranges  of  forces  and  moments, 
and  proper  balance  selections  may  be  made  for  accuracy  at  small  angles  of 
attack  or  for  a  range  to  cover  the  larger  angles  of  attack.  Model  support 
in  the  tunnel  is  positioned  so  that  dual  windows  can  provide  complete  flow 
visualization  and  photography.  Each  window  is  equipped  with  either  a  con^- 
ventional  single  pass  Schlieren  system  or  a  parallel  beam  refocused  shadow¬ 
graph  system.  Model  attitude  in  the  tunnel  was  measured  by  the  data  readout 
system  and  verified  optically  with  a  scope  in  both  the  pitch  and  yaw  planes. 

The  angle  of  attack  was  maintained  to  within  +  n.l<>  for  this  entire  inves¬ 
tigation. 

2.2  TEST  M0I«L 

The  model  used  in  this  investigations  was  a  7.25^  half  angle  cone  with  a 
tangent  full  dome  afterbody.  The  cone  base  diameter  is  10.6  inches.  A 
sketch  of  the  model  is  presented  in  Figure  2. Three  thermocouples  are  installed 
in  the  model  one  in  each  porous  chamber.  Additional  pressurd  taps  are  located 
inside  the  model  for  use  in  calibration  prior  to  the  test. 

The  model  is  longitudinally  divided  into  four  separate  sections.  The 
nose  section  is  fabricated  from  stainless  steel.  Sections  one  and  two  are 
fabricated  from  a  porous  sintered  nickel-stainless  steel  alloy  which  is 
approximately  60  percent  dense.  The  rear  plate  of  the  last  conical  section 
is  adapted  to  accommodate  the  hollow  compression-strut  support  system  used 
for  the  pressure  measurement  phase.  The  last  section  la  the  impervious 
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tangent  aft-dome.  It  is  made  of  stainless  steel  and  was  spin  fozned.  This 
section  has  cutouts  in  it  to  accommodate  the  legs  of  the  strut  support  system. 
A  second  dome  section  had  a  hole  in  the  center  to  accoiModate  the  sting  mount 
system,  which  was  used  for  force  measureinents. 

The  Internal  structure  of  the  model  is  fabricated  to  house  the  balance 
water  jacket  and  sting  arrangement  to  be  used  during  the  force  measurement 
phase.  The  internal  pltubing  of  the  model  is  designed  to  exhaust  gas  into  a 
coomon  plentsn  at  the  rear  of  the  model  inside  the  last  section.  The  gas  is 
dispersed  forward  through  sonic  orifices  into  the  two  porous  sections  of 
the  model.  The  mass  addition  rate  is  controlled  by  the  plentn  pressure 
and  the  distribution  of  the  mass  addition  is  controlled  by  the  relative  areas 
of  the  sonic  orifices  to  the  individual  chambers. 

The  blowing  or  porosity  uniformity  of  the  model  was  established  prior 
to  the  test  with  a  venturi  calibrating  device.  The  calibration  was  performed 
with  the  external  surface  of  the  model  at  atmospherid  conditions.  It  was 
performed  with  a  conical  venturi  which  had  a  twenty  to  one  contraction  ratio 
with  a  throat  area  of  0.05  square  inches.  The  differential  pressures  between 
throat  and  plentxn  were  recorded  with  a  +  1  PSIO  transducer.  The  gas  samples 
were  obtained  by  attaching  a  contoured  nibber  adaptor  which  was  fitted  to  the 
cone  surface.  The  capture  area  of  the  adaptor  was  roughly  0.5  square  inches. 
The  venturi  was  calibrated  with  the  AEDC  faciltiy  reference  sonic  orifices. 
Localized  mass  flux  readings,  (p  v)^,  were  obtained  along  several  rtyt  of  the 
model  which  were  45  degrees  apart.  The  calibration  indicated  that  the  model 
porosity  was  symmetric  to  within  t  25Z.  A  sample  of  the  calibration  data  is 
given  in  Figure  3. 

2.3  BOUNDARY  LAYER  PROBE 


A  pitot  pressure  and  a  total  temperature  probe  were  used  to  determine 
the  details  of  the  boundary  layer  during  Che  course  of  this  investigation. 

The  probing  apparatus  was  designed  to  present  both  prc^s  to  the  flow  with 
a  separation  distance  of  approximately  1/2  inch.  The  probing  was  performed 
with  each  probe  1/4  inch  off  of  Che  vertical  meridian  of  the  model  at  all 
angles  of  attack.  This  offset  was  accounted  for  in  the  data  tabulations  and 
plots.  The  probes  and  their  supporting  structure  were  incliided  7.25  degrees 
to  the  centerline  to  minimize  interference  with  the  flow  over  the  cone.  The 
temperature  probe  was  a  singly  shielded  thermocouple  with  a  diameter  of  0.060 
inch.  The  total  head  probe  was  a  .040  inch  diameter  hyperdermic  tube. 
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SECTION  III 

EFFECT  OF  BLOWING  RATE  NON-UNIFORMITIES  ON  BOUNDARY  LAYER  PROPERTIES 


I 


Prior  to  discussing  the  enalysis  of  the  data  obtained  with  the  porous 
cone  described  in  Section  ?,  some  words  are^in  order  with  regard  to  the 
sensitivity  of  experimental  data  to  the  nonuniformities  in  porosity  which 
are  indigeneoua  to  each  particular  model* 

Aerodynamic  test  models  fabricated  from  porous  materials  have  commonly 
been  used  to  simulate  ablation  in  wind  tunnel  experiments.  Diese  models  have 
been  constructed  from  a  variety  of  porous  materials  that  are  ctirrently  available 
commerically.  For  impulse  facilities*  such  as  AEDC  tunnel  F  (hotshot)  LTV* 
and  at  GE*  thin  skinned  modelf  (.030"  to  .060")  were  generally  fabricated 
(roiled  and  welded)  from  material  supplied  by  Union  Carbide,  Reendan  Associates* 
&Din  the  Pall  Corporation.  For  the  more  severe  heating  loads  and  environnent 
imposed  on  models  tested  in  the  continuous  flow  facilities  at  AEDC  tunnels  B 
and  C*  these  thin  skin  models  are  unsuitable.  In  these  facilities*  models 
fabricated  from  a  porous  material  frith  a  thicker  wall  (~.2S")  are  necessary. 
Materials  such  as  the  sintered  metal  of  Mott  Metallurgical  Corporation  or  the 
Foametal  material  of  GE  Metallurgical  Products  Division  are  therefore  required 
from  a  thermal  strese  viewpoint. 

Regardless  of  the  porous  material  selected*  the  porosity  distribution  - 
that  is  the  blowing  rate  distribution  over  the  model  -  must  be  known  for  each 
experiisental  investigation^  For  all  the  materials,  the  forming  process  in 
making  the  bulk  material  into  a  test  model  creates  certain  nonuniformities  in  the 
bloving  rate.  If  a  porous  model  Is  to  be  used  to  investigate  the  boundary 
layer  details  in  a  turbulent  mixing  investigation  for  example*  one  must  as¬ 
certain  what  nonuniformities  are  tolerable  from  a  data  accuracy  point  of  view. 

The  purpose  of  this  section  is  to  demonstrate  vith  experimental  data  and 
also  with  a  theoretical  calculation  conducted  vith  the  Equilibrium  Nonsimilar 
Boundary  l^i^er  Program  (ENSBL)'"^^,  the  insensitivity  of  the  mean  boundary 
Ij^^er  properties  to  the  ds-.'iationi;  from  the  everage  bloving  rate. 

3.1  NIMERKAL  VERIFICATION 


The  Equilibriun  Non-Similar  Boundary  Layer  Program  (ENSBL)  was  used  to 
investigate  the  effect  of  Che  distribution  of  mass  addition  race  on  the 
character  of  the  boundary  layer.  Calculations  were  carried  out  on  a  sharp 
7.25  degree  cone  for  Mach  nunber  8.  It  is  evident  (Figure  3)  that  the  blowing 
distribution  for  this  model  is  not  constant  but  varies  some  ±  50Z  about  the 
mean  value,  with  -f  25%  variations  in  axisymmetry. 

Ntmterical  results  were  cenputed  for  both  the  uniform  and  the  nonuniform 
(sawtooth)  axisymmetric  blowing  disttlbutions  shown  in  Figure  4.  In  one  case, 
a  uniform  mass  addition  rate  was  applied  over  the  porous  portion  of  Che  cone 
(0.19  ^  X/L^0.90).  0ut..ide  of  this  range  the  mass  addition  rate  was  zero. 
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I'n  a  second  case,  the  "sawtooth"  distribution  of  mass  addition  rate  was  used. 

This  distribucion  corresponds  to  a  fairing  of  the  data  of  Figure  3.  Total 
mass  added  was  the  same  in  the  two  cases.  Representative  results  are  plotted 
in  Figures  4  through  6. 

Qualitatively,  the  displacement  thickness  responds  as  expected  to  the 
mass  addition  rate.  Initially  (0.66  £  X  ^  1.18  ft.)  the  mass  addition  rate 
is  greater  in  the  uniform  ift  case,  as  is  the  displacement  thickness.  For  1.18 
ft.  ^  X  ^2.24  ft.  the  sawtooth  m  exceeds  the  uniform  value,  and  5*  for  the 
sawtooth  case  grows  faster  than  for  the  uniform  case.  At  X  ^  2.59  ft  and 
X  >  3.14  ft. the  mass  added  to  the  boundary  layer  forward  of  that  station  is 
the  same  in  the  tvo  cases,  as  is  the  value  of  displacement  thickness. 

Quantitatively,  the  response  of  the  curves  is  less  than  proportional 
to  the  variations  in  mass  addition.  At  X  1.9  ft.,  where  the  .mass  addition 
rate  for  the  sawtooth  case  is  537.  greater  than  that  of  the  uniform  case,  the 
difference  in  displacement  thickness  is  only  6Z.  At  X  =  2.2,  where  mass  added 
to  the  boundary  layer  in  the  sawtooth  case  is  15^  greater,  the  difference  is 
§  *  is  only  8%. 

Figures  5  and  6  present  velocity  and  stagnation  temperature  profiles  at 
X  =  2.2  ft. (where  the  maximum  difference  in  mass  added  occurs)  and  at  X 
3.065  ft.  (where  total  mass  added  in  the  two  cases  is  the  same).  As  might  be 
expected,  the  difference  between  the  two  cases  is  greater  in  Figure  5.  Rowver, 
the  differences  between  the  two  cases  is  insignificant*.  The  maxinun  difference 
in  temperature-  profiles  is  less  than  17,,  whereas  the  maxinun  difference  in 
velocity  profiles  is  less  than  27,. 

The  conclusions  ^ich  can  be  dratm  from  this  comparison  are;  because 
of  the  increased  mixing  characteristics  in  a  turbulent  boundary  layer,  a 
difference  in  the  mass  addition  rate  distribution  about  a  mean  value  does 
not  appear  to  significantly  affect  the  character  of  the  boundary  layer;  and 
that  the  boundary  layer  depends  primarily  on  average  mass  addition  rate,  and 
appears  to  be  relatively  insensitive  to  fluctuations  in  the  upstream  history 
of  the  mass  injection. 

3.2  EXPERIMENTAL  VERIFICATION 

Force  and  moment  data  were  obtained  with  the  porous  model  mounted  on  a 
strain  gage  balance.  Data  were  obtained  for  several  blowing  races.  If  the 
blowing  nonuniformities  cause  an  asymmetric  configuration  to  result,  then  a 
residual  moment  and  normal  force  would  appear  at  ot  =  0°.  However,  data 
obtained  with  this  model,  as  shown  in  Figure  7,  show  no  such  influence. 

It  car.  be  concluded  that  useful  and  meaningful  experiments  can  be  con> 
ducted  to  investigate  the  turbulent  boundary  layer  characteristics  of  an 
ablating  cone,  simulated  by  blowing,  with  a  model  which  exhibits  blowing 
nonuniformities  of  +  507,  about  the  mean. 


It  should  be  nrted  that  a  highly  expanded  temperature  scale  was 
used,  i.e.,  only  the  range  0.8^Tq/Tq  fil.O  is  shown. 

OB 
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SECTION  IV 


ANALYSIS  AND  CORRELATION  OF  TURBULENT  BOUNDARY  LAYER  DATA 


The  experimental  results  obtained  from  the  boundary  layer  profile  measure¬ 
ments,  which  include  the  effects  of  mass  injection,  are  presented  and  discussed 
in  the  following  sections.  The  data  analyzed  are  for  both  turbulent  and 
laminar  boundary  layer  flows,  although  emphasis  is  placed  on  the  turbulent  data. 
For  the  turbulent  data,  the  effect  of  mass  injection  on  such  properties  a.~ 
the  boundary  layer  thickness,  displacement  and  momentum  thicknesses,  velocity 
profile  exponent  and  the  applicability  of  the  Crocco  theory  are  discussed. 

The  data  analyzed  are  presented  in  graphical  form  in  the  body  of  the  report; 
tabulations  of  the  profile  properties  deduced  from  the  measured  data  are 
contained  in  the  Appendix  A.  Also  included  in  this  report  are  the  free  stream 
and  local  properties  ti*^*  wall  temperature,  pressure,  transition  location, 
etc.)  that  are  associated  with  the  profile  measurements. 

4.1  DATA  ANALYSIS  AND  PRESENTATION 

The  boundary  layer  profile  data  that  were  obtained  with  the  porous  cone 
shown  in  Fig.  2  for  both  the  1969  STREET-G  program^^^)  and  also  for  the  1968 
Mark  12  Penetration  Aids  program(28)  win  be  presented  in  this  section.  For 
convenience  and  clarity  of  presentation,  the  data  were  coded  to  denote  the 
free  stream  environment,  bloteing  rate  and  distribution,  the  injectant,  and 
the  state  of  the  boundary  l^er  at  the  profile  measurement  station.  Each 
case  is  assigned  a  letter  index  lAich  denotes  the  free  stream  condition, 
the  average  blowing  rate  and  distribution,  and  the  state  of  the  boundary  layer 
at  the  measured  profile  station  (see  Table  I).  A  single  letter  denotes  the 
data  acquired  on  the  1969  STREET  G  Program'^”).  A  double  letter  denotes  data 
acquired  cn  the  1968  Mark  12  ?/A  effort(^);  moreover,  the  second  letter  L 
denotes  laminar  data.  Each  case  is  assigned  a  data  symbol  which  will  be  used 
consistently  throughout  the  report  (see  Table  I) .  For  each  case  surface 
pressure,  force  data,  and  boundary  layer  profile  data  at  several  stations 
were  obtained  from  Reference  26.  The  AEDC  data  group  numbers  corresponding 
to  each  case  are  shown  in  Table  2.  Average  free  stream  test  conditions  for 
each  case  may  be  found  In  Table  3. 

4.1.1  WALL  TEMPERATURE  DISTRIBUTION 

The  model  wall  temperature  distributions  were  determined  for  each  blowing 
condition,  at  each  Reynolds  nunber.  The  temperature  distributions  as  determined 
from  the  three  thermocouples  located  in  the  porous  sections  are  shown  in 
Figure  8.  These  values  correspond  to  the  equilibriun  wall  temperature  that 
is  achieved  prior  to  and  maintained  during  the  acquisition  of  the  profile  data. 
For  each  Reynolds  nunber,  the  wall  temperature  as  deduced  from  each  thermo¬ 
couple  correlates  with  the  local  average  blowing  rate,  and  as  sho.’i  in  Figure  8 
is  reduced  as  the  blowing  rate  increases.  This  was  also  true  fr'r  the  non- 
uniform  blowing  cases  (M  and  N) . 
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4.1.2  TRANSITION  ZONE  DEFINITION 

Measurements  to  deduce  the  state  of  the  boundary  layer  were  not  made  during 
the  1969  STREET-G  tests,  primarily  because  these  measurements  were  available 
from -the  prior  test  entry The  state  of  the  boundary  layer  was  determined 
from  the  impact  pressures  measured  by  traversing  a  pitot  probe  along  the  model 
surface.  Data  to  assess  the  effect  of  mass  addition  on  boundary  layer  transi¬ 
tion  were  obtained  for  one  test  condition  (i.e.  Re  „  /ft  =  2.7  x  10°)  and  are 
shown  in  Figure  9.  Variation  of  onset  and  end  of  transition  with  mass  addition 
rate  were  deduced  from  Figure  9,  and  are  plotted  (solid  lines)  in  Figure  10. 

The  location  of  the  transition  zone  for  the  second  test  condition  (i.e,  Re^/pt.  = 
3.7  X  10  )  was  obtained  for  the  non  blowing  case  only.  The  m  effects  at  the 
higher  Reynolds  number  cc'..dltlon  (dotted  lines)  were  then  scaled  from  thii  data 
available  at  the  lowc.  ;^:ynolds  number.  It  is  apparent  from  the  dip  (near 
m  =  .01)  in  the  data  shown  in  Figure  10,  that  a  slight  transition  reversal 

occurred.  The  initial  forward  movement  can  be  attributed  to  the  destabilizing 
influence  of  mass  addition  on  transition.  However,  as  the  blowing  rate  Increases, 
the  wall  temperature  decreases.  In  this  temperature  range,  the  litter  effect 
is  stabi.  .zing,  and  apparently  is  dominant,  so  that  as  the  blowing  rate  in¬ 
creases  beyond  0.01  transition  moves  aft.  Furthermore,  the  data  Indicate  that 
the  local  wetted  length  Reynolds  nianber  (Res)  at  the  end  of  transition  is 
roughly  8-9  x  10”.  In  a  later  discussion  concerning  the  data  of  Danberg^^^) 
the  above  information  will  be  used  to  indicate  that  some  of  the  data  presented 
therein  are  not  turbulent  but  are  transitional. 

4.1.3  MASS  ADDITION  RATE  DISTRIBUTION 

Consistent  with  the  calibrations  performed  with  the  porous  model  at  AEDC, 
the  mass  addition  rate  distributions  were  scaled  to  confarm  with  the  basic 
calibration  and  ciie  integrated  (metered)  mass  rate  to  the  model.  Listed  iu 
Table  4  are  the  mass  addition  rate  schedules  for  Cases  A  through  N. 

4.1.4  BOUNDARY  LAYER  PROFILE  MEASUREMENTS 

As  i .  -icated  in  Section  2.3,  boundary  layer  profile  measurements  (i.e. 
pitot  pressure  and  total  temperature)  were  made  at  several  axial  stations 
(see  Figure  2).  The  Rayleigh  formula^^®)  was  used  to  compute  the  Mach  number 
from  the  measured  pitot  pressure  and  the  static  pressure  which  was  assumed 
constant  in  the  boundary  layer  and  equal  to  the  wall  static  pressure,  llie 
wall  static  pressure  was  measured  by  orifices  in  the  model  surface  in  the 
absence  of  the  profile  probes.  The  velocity  in  the  boundary  layer  was  com¬ 
puted  from  the  local  value  of  the  Mach  number  and  the  measured  total  tempera¬ 
ture.  Using  adiabatic  relations^^®)  in  this  manner  boundary  layer  properties 
such  as  density,  temperature,  and  velocity  which  are  required  to  compute  the 
displacement  and  momentum  thicknesses  were  defined.  The  displacement  and 
momentum  thicknesses  were  deteimlned  from  the  following  relations; 
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Inherent  in  the  above  relations  are  the  corrections  which  account  for 
transverse  curvature.  The  boundary  layer  thickness  was  established  consistently 
from  the  total  temperature  profile,  which  for  the  sharp  cone  also  agrees  with 
that  obtained  from  the  pitot  pressure  profile.  Tabulations  of  the  local  con¬ 
ditions  and  profile  parameters  may  be  found  in  Appendix  A.  For  Re^o  /ft  = 

3.7  xlO°  (cases  A  through  F.)  the  axial  distribution  of  6,  5  *,  and  0  may 
be  found  in  Figures  11  through  13,  and  for  Re^/ft  =  2.7  x  10®  (cases  F  through 
J)  in  Figures  14  through  16.  The  increase:  in  the  viscous  thickness  due  to 
mass  injection  is  rather  evident  from  these  figures.  It  will  be  shown  that 
mass  transfer  Increases  each  thickness  roughly  linearly  with  the  blowing  para¬ 
meter  (^v)  *  Also  shown  in  these  figures  is  the  boundary  layer 

j  -  ^ 


u), 


w 


transition  zone  for  m  »  0.  As  noted,  the  forward  profile  station  (x/L  =  0.35) 
is  in  the  transition  zone.  Also  shown  on  these  figures  is  the  zero  mass  addi¬ 
tion,  all  turbulent  boundary  layer  prediction  of  WalkerC^^),  as  computed  by 
V1ZAAD'^^\  For  the  sharp  cone.  Walker’s  relations  are  simply  defined  as 


e 


0.  037  (MF)  - 


(3) 


vdiere  MF  is  the  Mangier  factor 


<4) 

(*) 


where  equations  4  and  5  are  valid  for  Me  2  4.5.  It  is  evident  that  modifi¬ 
cations  must  be  made  to  these  simple  engineering  relations  for  the  zero  mass 
addition  case  in  order  to  account  for  the  influence  of  mass  injection.  This 
will  be  discussed  in  Section  5.1. 

Turbulent  boundary  layer  velocity  profiles  are  cv  '.tomarily  characterized 
by  the  power-law  relation 

(•i)  " 
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The  exponent,  nj,  was  obtained  by  measuring  the  slope  of  the  velocity  profile 
plotted  on  log-log  paper.  The  axial  distribution  of  n  for  Cases  A  to  E  is 
shown  in  Figure  i7,  and  for  Cases  F  through  J  in  Figure  18.  For  the  zero 
injection  case  at  /ft  =  3.7  x  10°,  the  axial  distribution  is  characterized 

by  a  value  of  n  »  9  at  the  forward  transitional  station,  an  increase  to  sene 
unknown  value  (.n  12)  at  the  end  of  transition,  then  by  a  decrease  to  values 

of  the  order  of  7  at  points  sufficiently  far  removed  from  the  end  of  transition. 
This  is  in  concert  with  the  study  performed  by  Johnson  and  Bushnell^^^)  wherein 
a  data  summary  concerning  the  power  law  sensitivity  to  Reynolds  number,  Mach 
number,  and  wall  cooling  was  made  (Figure  19).  The  Johnson-Bushnell'^*'  data 
survey  indicated  that  the  velocity  exponent  actually  experiences  an  overshoot 
beyond  the  classical  turbulent  value  of  n  «»  7  at  points  immediately  downstream 
of  transition.  It  was  postulated  that  n  correlates  with  the  incremental  dis¬ 
tance  downstream  of  the  end  of  transition  normalized  by  the  boundary  layer 
thickness  (/i.X^^/5).  However,  the  present  data  do  not  appear  to  support 
this  hypothesis.  Nonetheless,  values  of  the  velocity  exponent  in  excess  of 
seven  were  experienced  in  the  region  immediately  downstream  of  transition. 


The  effect  of  mass  transfer  on  the  turbulent  velocity  profile  shape  is 
to  make  it  less  full  and  as  a  result  the  exponent  n  is  reduced  as  shown  in 
Figures  17  and  18.  A  decrease  of  n  to  some  30%  of  its  no  blowing  value  is 
noted  for  a  blowing  rate  C  ^  • 

4.2  DATA  CORRELATIONS  AMD  COMPARISONS 

To  assess  the  Influence  of  mass  addition  on  the  viscous  layer  thicknesses, 
the  turbulent  data  available  in  Refs.  26  and  28  were  presented  in  a  form  where 
the  normalized  thickness  ratio  (blowing  value  to  the  no  blowing  value)  were 
plotted  against  the  blowing  parameter,  The  values  of  Cfo  and  (P  u)eo 

represent  local  properties  at  the  profile  station  for  the  non-blowing  case. 


The  thickness  growth  with  blowing  for  6/6 


6*/ 6*  and  0/6  are  showh 
o  o 


in  Figs.  20  to  22,  reflectively.  Also  included  in  these  figures  are  the  tur¬ 
bulent  flat  plate  data  of  Danberg(29)  where  the  ^  value  was  corrected  by  the 
Mangier  factor.  Good  agreement  and  correlation  is  noted  for  all  the  data. 

As  mentioned  earlier,  the  data  presented  by  Danberg  correspond  to  data  obtained 
for  two  separate  ranges  of  wetted  length  Reynolds  number,  Reg>  one  covering  the 
range  Res  of  3  x  10°  to  6  x  10^,  and  the  second  covering  the  range  8.5  x  10® 
to  11  X  10®.  It  is  postulated,  that  based  upon  the  available  boundary  layer 
transition  results,  that  the  data  obtained  for  the  lower  Reynolds  number  range 
correspond  to  transitional  boundary  layer  data.  When  these  data  are  included 
in  the  correlation,  although  agreement  is  achieved  with  the  mean  levels  of 
the  "turbulent"  data,  considerably  more  scatter  is  evident.  Consequently,  in 
the  correlations  shown  in  Figs.  20  to  22  only  data  for  P^es  ^  9  x  10®  are 
included.  Best  fit  curves  were  established  for  each  of  the  ratio's  in  the 
form  of  simple  algebraic  relations.  That  is*. 

6/6^  “  1  +  0,2(52  C-  0.00627 

d*/6*  »  1  +  0.333? 
o 

e/Sp  -  1  +  0.410  c 
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As  a  result  of  these  curve  fits,  the  effect  of  mass  addition  on  the  form 
factor  H  1  6  */6  becomes 

H  »  I-*- 0.333  C 
H-  1  + 0.410  C 


A  comparison  of  this  relation  and  the  data  is  shown  in  Fig.  23.  It  is  evident 
from  these  data  that  mass  addition  tends  to  lower  the  form  factor  so  that  at 
C  ^  9,  there  is  roughly  a  15%  This  is  in  contrast  to  all  of  the 

existing  turbulent  correlations^'^*^^*^^'  which  assume  that  blowing  has  no 
effect  on  the  form  factor  (i.e.  H  =  Hq). 


Also  shown  in  Figs.  20  to  22  are  the  curves  corresponding  to  the  current 
empirical  prediction  techniques  of  General  Electtic^'^)  McDonneliDouglas(32), 
and  Phllco  Fovd(21).  jhe  analytic  form  for  each  of  these  relations  is  ex¬ 
pressed  below: 

General  Electric^ 


.  1  r  4  +  4T  e 

“T  TTT 

L  J  o 


6V6* 

6*/6* 

o 


1  +  0.752  C  +  0.0189  C' 
0. 185  +  t 


McDonne ll Douglas^ ^ 


(0  s:  C  s  6.55)  (12) 
(C  >  6.55) 


1  +  0.143^ 


W®*  ®  C 

6V6*  -  6/0  »  - 0  8'  "? — 


^erew=  Vin/rb.l. 


is  the  ratio  of  specific  heats 


Phllco  Ford^^^) 


6*/6* 

o 


•  1  + 


-- (-f-)'”"  (tr 

1/2 

1  /qr,0.46  -  ^\2  1  /~  0.46  ~  ~  ,.\2  1 

Ti“w  "  «)  -t(“w  “')[(“w 
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It  Is  evident  from  the  data  comparisons  that  each  of  the  techniques 
described  above  overpredicts  the  effect  of  mass  addition  on  the  viscous  layer 
thicknesses. 

The  effect  of  mass  addition  on  the  velocity  profile  exponent  is  shown 
in  Fig.  24,  It  can  be  seen  from  this  figure  that  a  considerable  distortion 
of  the  velocity  profile  results  from  mass  addition.  Furthermore,  the  scatter 
of  the  results  suggests  that  parameters  in  addition  to  ^  are  necessary  to 
better  correlate  the  data. 

4.2.1  SENSITIVITY  OF  RESULTS  TO  AJ^GLE  OF  ATTACK 

It  is  desirable  to  assess  the  sensitivity  of  the  data  used  in  the  cor¬ 
relations  to  various  sources  of  error,  for  example  to  the  model  angle  of 
attack.  This  is  particularly. true  of  the  results  presented  in  Figures  20 
through  24,  since  experimental  error  is  magnified  in  obtaining  ratios  of  two 
measured  quantities.  Although  deliberate  care  was  taken  to  align  the  model 
precisely  with  respect  to  the  oncoming  stream  (i.e.Acx  *  +0.1°)  some  mis¬ 
alignment  is  always  present.  Boundary  layer  profiles  were  obtained  on  ghe 
windward  and  leeward  rays  of  the  model  at  selected  stations  for  a  =fl 
(see  Ref.  26).  A  data  reduction  procedure  similar  to  that  used  for  the 
a  =  0°  cases  was  used  with  the  C:  =  +  1°  zero  mass  addition  data  at  x  = 
41.51  inches  (chat  is, cross  flow  effects  were  neglected).  Shown  in  Fig.  25 
are  the  results  of  this  calculation.  An  error  of  +  2X  in  thickness  would 
result  for  a  ^  0.1°  angle  of  attack.  However  the  effect  of  angle  of  attack 
on  the  velocity  profile  exponent  is  negligible,  to  within  the  accuracy  of 
the  data.  Consequently  the  data  scatter,  evident  in  Figure  24,  must  be 
attributed  to  other  causes. 

4.2.2  COMPARISONS  OF  D.ATA  WITH  THE  CROCCO  RELATION 

Theoretical  and  semi-empirical  approaches  to  the  solution  of  a  turbulent 

boundary  layer  rely  upon  several  approximations,  one  of  which  is  the  Crocco 

relation  for  unit  Prandtl  number,  H  -  h  The  bulk  of  available 

w  _ 


experimental  data  indicate  that  only  for  the  adiabatic  wall  case  does  the 
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unit-Prandtl  number  Crocco  relation  agree  with  the  data.  In  terms  of  the 
total  enthalpy  profile,  the  sunmary  of  experimental  results  reported  by 
Bertram  and  Neal^^^)  indicates  consistently  lower  values  of  total  enthalpy 
in  the  turbulent  boundary  layer  than  the  Crocco  relation  (see  Fig.  26). 

Purely  as  an  empiricism,  a  curve  for  H  -  h  /  u  \  ^  also  drawn  in 

H  -h  “  ('ir) 

e  w  \  e  / 

Fig.  26, and  is  seen  to  describe  much  more  adequately  the  experimental  results 
than  the  classical  linear  Crocco  relation.  Included  in  this  figure  are  the 
results  of  the  present  investigation^^^) .  For  the  present  data,  the  absence 
of  data  at  u/u^  <  .7  ie  due  to  the  difficulty  of  obtaining  accurate  measure¬ 
ments  In  the  relatively  thin  boundary  layers  present  on  slender  cones  in 
hypersonic  flows.  Nonetheless,  the  present  data  for  the  non-blowing  case 
agree  with  the  data  summary  of  Bertram  and  Neal(33). 


The  effect  of  blowing  on  the  total  enthalpy  profile  for  all  stations  in 
case  E  may  be  seen  in  Fig.  27.  In  this  case,  due  to  the  boundary  layer 
thickening  resulting  from  mass  addition,  data  were  obtained  down  to  u/ug  = 
0.27.  Furthermore,  the  accuracy  of  the  measurement  (reduction  of  scatter) 
is  also  shown  to  increase.  It  is  evident  from  this  figure  that  for  a  given 
value  of  u/ug  mass  injection  tends  to  increase  the  total  enthalpy  in  the 
boundary  layer  over  the  zero  m  value.  A  consistent  trend  was  noted  for  all 
of  the  data  obtained,  vdierein,  an  increase  in  blowing  affected  the  level  of 
the  total  enthalpy  in  the  boundary  layer  (Fig.  2b),  Thus,  if  one  represents 
the  enthalpy  profile  by  the  relation 

(18) 


H  -  h 
H 


•  h 

-h  \  u 
w  '  e 


the  value  of  €  (shown  in  Fig.  29)  tends  to  decay  with  increased  mass  addition. 
It  was  tacitly  assumed  that  €=  2  best  represents  the  zero  mass  addition  data. 

4.3  DATA  COMPARISONS  WITH  THEORY 


Comparisons  of  the  boundary  layer  thickness  data  with  the  General  Electric 
Theories,  such  as  the  Viscous  Interaction  Zero  Angle  of  Attack  Drag  (VIZAAD), 
and  the  Equilibrium  Non  Similar  Boundary  Layer  (SNSBL)  Programs,  will  be  made 
in  this  section.  The  techniques  employed  in  the  current  versions  of  each  of 
these  programs (^®> 27,20)  ^j-g  nominally  equivalent  to  those  employed  industry 
wide  and,  consequently,  the  comparisons  serve  to  show  where  deficiencies  in 
the  current  state  of  the  art  exist.  In  addition,  the  implication  of  the 
recommended  changes  to  the  VIZAAD  program  on  flight  data  predictions  will  be 
assessed. 

4.3.1  CCMPARISONS  OF  DATA  WITH  VIZAAD  PROGRAM 

Predictions  of  the  viscous  layer  thickness  as  generated  by  the  VIZAAD 
program  were  compared  with  data  cases  A  through  E  as  shown  in  Figs.  30  to  34. 
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In  these  VIZAAD  computations,  the  flow  was  treated  as  laminar  up  to  the  point 
of  transition  onset.  Downstream  of  this  point,  the  flow  was  treated  as  tur¬ 
bulent.  The  effective  origin  of  the  turbulent  solution  was  established  by 
matching  the  zero  blowing  laminar  momenttmi  thickness  to  the  zero  blowing  tur¬ 
bulent  momentum  thickness,  that  is  (0)  lAMXNAR  *  TURBULENT*  ^  con¬ 
sequence  the  values  of  the  boundary  layer  tnickness,  5  ,  and  the  displacement 
thickness,  5  are  discontinuous  at  the  transition  point  even  for  the  zero 
mass  addition  case.  For  the  cases  with  mass  injection,  all  three  thicknesses 
suffer  the  discontinuity  since  the  thickness  correlations  differ  at  this  poii.t. 
For  the  cases  under  consideration,  transition  onset  was  sufficiently  far  for¬ 
ward  so  that  the  influence  of  this  discontinuity  on  the  downstream  values  is 
small.  It  Is  evident  from  Figure  30  that  the  original  correlations  grossly 
underpredict  the  viscous  layer  thicknesses.  Implicit  In  the  constant  of 
equation  3  for  computing  the  momentum  thickness.  Is  the  assunption  that  the 
velocity  profile  exponent  n  in  the  entire  turbulent  regime  18  7.  For  a  flat 
plate,  the  origin  of  Ll;“  '  stant  may  be  found  in  Schlichting^^)  which  can 
be  summarized  as  follows: 


For  the  general  case  of  the  power  law  u/u^ 


(y/6)'^" 


we  have 


J.. 

n 

6 

(1  +n)(2+n) 

and 

6  = 

/  \®*^ 

0.0575 

X 

J  (l  +  n)(2+n)  j 

0  ** 

(19) 


0.8 


(20) 


which  for  n 

=  7 

becomes: 

,  .-0.2 

»  0. 37  X  ( 

(21) 

6 

<  V  1 

and 

/  \  —0. 2 

0 

»  0.036X  1 

(22) 

^  1/  ) 

However,  the  experimental  value  of  n  is  larger  than  the  theoretical  value 
predicted  by  Equation  3.  This  implies  that  n  should  be  greater  than  7,  which 
is  what  the  values  deduced  from  the  profiles  indicate.  If  a  functional  rela¬ 
tion  n  (x)  were  known,  the  constant  could  be  revised  accordingly,  through  the 
use  of  equation  20.  However,  this  is  not  generally  known.  From  the  data 
of  Cases  A  and  F  (that  is  for  m  =  0)  one  ran  deduce  that  the  constant  In  the 
momentum  thickness  relation  should  be  increased  from  0.037  to  0.059.’  The 
displacement  and  the  boundary  layer  thickness  will  increase  proportionately 
as  is  evident  from  equations  4  and  5.  The  revised  relations  show  good  agree¬ 
ment  with  the  data  as  seen  in  Fig.  30. 
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since  the  effects  of  mass  addition  on  the  thickness  are  currently  being 
overpredicted,  the  existing  program  tends  to  over  compensate  for  the  poor 
agreement  at  the  lower  blowing  rates  as  is  evidenced  by  Figs.  31  to  34.  The 
relations  defining  the  effects  of  mass  addition  in  the  program  (that  is 
equations  11  and  12)  were  replaced  by  the  modified  relations  (equations  7 
through  9).  Shown  in  Figs.  31  to  34  are  comparisons  of  the  revised  program 
with  the  data.  It  is  noted  in  ail  cases,  that  good  agreement  was  achieved  by 
all  three  thicknesses. 

For  blunt  bodies,  the  properties  at  the  edge  of  the  boundary  layer  are 
determined  (in  VIZAAD)  from  the  mass  balance  relation 


m  *  m. .  +  m 

CO  inj.  b.  1. 


(^u)  ir  •  If  (^*v)  dA+2ff  r.  (P\i)  (6-6*) 
<*>  s  J  J  w  b  e 


It  is  evident  that  gross  modifications  to  the  thickness  values  will  affect 
this  balance  and  accordingly  will  affect  the  local  properties  and  conse¬ 
quently  vehicle  drag.  Flight  data  comparisons  for  two  classes  of  slender 
vehicles  were  made  and  some  of  the  results  are  sumnarized  below.  Vehicle  1, 
is  a  small-mass-loss  ablator  (i.e*  has  a  high  temperature  heat  shield). 

Vehicle  2,  is  a  large-mass-loss  ablator  (i.e.  has  a  low  temperature  heat  shield), 
The  initial  bluntness  ratio  for  both  these  vehicles  isR^j, /Rg  =  .04. 


Comparison  of  Revised  Program  %rlth  Flight  Data 


Vehicle 

Altitude 

Vehicle 

Original 

Revised 

Flight 

No. 

Kft 

Station 

Parameter 

Version  of 

Version  of 

Data 

VIZAAD 

VIZAAD 

1 

25 

X  =  .5L 

6.  (ft) 

.0264 

.0288 

High 

6  (ft) 

.00813 

.00958 

Temp. 

Me 

3.50 

3.97 

Ablator 

Reg  X  10“6 

12.66 

16.57 

m  *  .001 

q  Ctu/ff^sec 

273.9 

283.2 

X  =  L 

5  *  (ft) 

.0169 

.0223 

6  *  (ft) 

.00746 

.0101 

5.67 

6.73 

Re  X  10~® 

124.5 

210.6 

4  Stu/ft^sec 

373. 

382.7 

— 

1 

Total  Vehicle 

.0724 

.0743 

.082 

Drag 

_ 

2 

80 

X  =  .5L 

6  .(ft) 

.0124 

.0164 

Low 

a  (ft) 

.00753 

,00857 

Temp. 

Me 

13.78 

14.30 

Ablator 
fe  =  .049 

Re.  X  10*® 
q  Btvi/  f  t* 

17.7 

598 

20.2 

579 

sec 

_ 

X  »  L 

a  (ft) 

.0431 

.0487 

a*  (ft) 

.0320 

.0287 

Me 

14.5 

14.14 

Reg  X  10*® 

4  Btu/ft^sec 

41.5 

38.1 

195 

189 

Total  Vehicle 

.0600 

.0570 

.043 

Drag 

For  these  comparisons  the  local  boundary  layer  edge  properties  such  as 
the  Mach  number  and  Reynolds  nunber  show  significant  variations  where  the 
bluntness  effects  are  felt.  Furthermore,  in  all  comparisons  with  flight  data, 
che  revised  version  of  the  program  predicted  drag  values  closer  to  the  data 
Irrespective  of  whether  the  predictions  '■«re  lower  than  or  higher  than  the 
flight  values. 
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Independent  of  the  drag  comparisons,  large  changes  in  the  viscous  layer 
predictions  result  from  the  modification.  It  will  be  shown  In  the  nej.t  section 
that  the  results  of  the  revised  VIZAAD  program  agree  batter  with  the  more 
exact  ENSBL  program  predictions. 

In  addition  to,  and  more  important  than  the  effect  of  the  thickness 
perturbation  on  local  properties,  is  the  effect  this  change  has  on  plasma 
attenuation.  Here,  there  is  interest  in  transmitting  at  higher  frequencies, 
tdiich  implies  shorter  wave  lengths,  lliese  wave  lengths  are  of  the  c'  der  of 
the  boundary  layer  thickness  and,  consequently,  the  thickness  of  tb'  layer 
and  the  shape  of  the  profile  (l.e.  the  local  value  of  n)  have  a  gros  '•  effect 
on  the  plasma  behavior.  Shown  in  Figure  35  is  the  S  band  plasma  attenuation 
profile  with  altitude  for  a  slightly  blunted  slender  re-entry  vehicle.  Also 
shown  on  this  figure,  at  altitude  h,  is  the  plasma  attentuation  as  computed 
by  the  method  of  Ref.  35  for  two  values  of  the  boundary  layer  thickness.  One 
was  deduced  from  Walker's  original  relations^^^)  and  the  second  was  computed 
from  the  relations  as  modified  in  this  report.  The  profile  shape  was  similar 
for  the  two  computations.  One  will  note  that  the  computation  is  sensitive 
to  the  thickness  of  the  boundary  layer,  and  furthermore  the  comparison  with 
data  is  better  with  the  modified  thickness  relations. 

4.3.2  COMPARISCWS  OF  DATA  WITH  THE  ENSEL  PROGRAM 

_  Comparisons  of  the  ENSBL  program  were  made  with  the  data  of  Case  A 
(A  =  0)  and  Case  C  (m  »  .01) .  In  the  former  comparison,  three  different 
eddy  viscosity  models  were  employed  in  a  manner  similar  to  that  discussed 
in  Reference  10.  For  the  latter,  comparisons  of  the  theory  with  one  eddy 
viscosity  model  were  made  with  the  measured  thicknesses. 

Shown  in  Fig.  36  are  the  results  of  the  ENSBL  calculation  compared  to 
the  data  of  Case  A.  The  three  eddy  viscosity  models  emgloyed  are  those 
proposed  by  Smith  and  Cebeci^^),  Patankar  and  SpaldingC®)  and  by  Sontowski^^)  . 

In  general  the  agreement  of  all  three  calculations  with  the  measured  total 
and  displacement  thicknesses  is  reasonable.  However,  one  will  note  that  the 
moroentum  thickness  is  underpredicted  by  the  program  irrespective  of  the  viscosity 
model  used.  Examination  of  the  velocity  profile  prediction  in  the  physical 
coordinates  (Fig.  37),  and  in  the  Crocco  coordinates  (Fig.  38)  indicates  that 
the  total  enthalpy  in  the  viscous  layer  is  being  overpredicted.  Of  the  three 
viscosity  models,  that  proposed  by  Sontovski  (which  employs  the  lowest  value 
of  turbulent  prandtl  nianber)  most  closely  resembles  the  data  (Fig.  38).  However, 
additional  modifications  are  required  to  effect  acceptable  agreement.  The  impli¬ 
cations  are  that  alterations  to  the  turbulent  Prandtl  number,  among  other 
things,  are  required. 

Comparisons  of  the  ENSBL  progra..,,  using  Sontowski's  eddy  viscosity  model, 
were  made  with  data  representing  a  case  with  mass  injection  as  shown  in  Fig. 

39.  Once  again,  reasonable  agreement  is  noted  for  the  total  and  displacement 
thickness  distributions.  However,  the  momentum  thickness  is  underpredicted 
as  in  the  zero  mass  injection  case.  One  can  speculate  that  although  mass 
injection  increases  the  enthalpy  in  the  viscous  layer,  the  theory  predicts 
larger  values  than  those  determined  in  the  experiment. 


17 


SECTION  V 


ANALYSIS  AND  CORRELATION  OF  LAMINAR  BOUNDARY  LAYER  DATA 

The  turbulent  data  reported  in  references  (26)  and  (28)  were  discussed 
in  the  preceding  section.  Nine  cases  in  which  the  boundary  l^er  flow  re¬ 
mained  laminar  (see  Table  1)  are  also  reported  in  Reference  (28).  The  purpose 
of  the  present  section  is  to  compare  the  laminar  data  obtained  in  these  cases, 
with  currently  used  empirical  correlations. 

5.1  DATA  ANALYSIS  AND  PRESENTATION 

The  boundary  layer  data  in  the  laminar  cases  were  obtained  and  reduced 
by  the  same  methods  as  described  in  Section  4.1  for  the  turbulent  cases. 
Displacement  and  momentum  thicknesses  were  computed  using  equations  (1)  and 
(2).  For  this  laminar  data,  tney,  as  well  as  boundary  layer  thickness  ,  6  , 
are  plotted  versus  axial  position  in  Figures  40  through  42. 

To  show  their  dependence  am  mass  addition  rate,  the  ratios  fi  /6  q, 

Q  /  and  a  I  f^Q  are  plotted  vs  C  (Figures  )3  th2:t>*jgh  45) .  Although 

cases  QL  and  SL  provided  data  for  values  of  C  near  30.0,  too  few  points 
were  available  at  these  high  mass  addition  rates  to  warrant  presentation. 
Figures  40  through  45,  were,  therefore,  limited  to  A.S  might  be  ex¬ 

pected,  the  data  show  a  general  thickening  of  the  boundary  layer  with  an 
increase  in  mass  addition  rate. 

5.2  DATA  COMPARISONS  WITH  CURRENT  THEORY 

Laminar  thickness  predictions  of  Walker^^^^  as  computed  in  VIZAAD^^^ 
are  compared  with  experimental  values  in  Figures  40  through  42.  Walker's 
relations  for  laminar  flow  are; 


S.  3  s 


u  “/2h  +3.36  h  /h  -  0.376 
0  *  _  e  8 _ w  e _ 

^  u  ^/2h  +  3.36  h  /b  +  5.79 
e  e  we 


and  Q  is  obtained  from  an  integration  of  the  following  equation,  which  is 
the  result  of  combining  the  integral  momentw  and  continuity  equations: 


d  2  a  f  X  6* 

ip  u  0  +  X  - - 

ds  e  2  u 

e 


du 

^  ^  _i _  « 

u  ds  2  ds 

e  0  u 

e  e 
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The  value  of  C  ,  which  is  required  in  equation  (26),  is  obtained  (in  the 
absence  of  mass ^addition)  through  the  Reynolds  analogy  from  Lee's  relations 
for  heat  transfer  rate,  modified  by  th'*  reference  enthalpy  method^^C). 


0. 354  P*  ti*  u  r.  (h  -h  ) 
8  b  r  w 


r  .s  1 

1/2 

p2/3 

1  P*M*u  r^dS 

r 

J  e  b 

L  o  J 

(27) 


At  the  higher  of  rhe  two  Reynolds  nunbers  for  wnich  laminar  data  were 
available  (Re*  =  2.7  x  10^/ft),  equations  (24)  and  (25)  consistently  over¬ 
predict  fi  and  5*  (Figure  40).  At  the  Icwer  Reynolds  number  (0.67  x  10^/ft) 
no  consistent  trend  can  be  observed.  Nor  can  one  deduce  any  trend  vdien  com¬ 
paring  the  momentum  thickness  given  by  equation  (26)  with  values  computed 
from  measured  data. 

Because  of  the  limited  amount  of  experimental  data  available,  no  attempt 
is  made  to  suggest  improvements  to  the  laminar  correlations. 


The  effect  of  mass  addition  on  momenttin  thickness  is  ccn.puted  in  VIZAAD 
by  first  solving  equation  (26)  for  ths  zero  mass  addif'on  case;  then  repeating 
the  solution  for  ^  ,  inserting  the  appropriate  of  (Pv)^  ir.  equation  (26)  ^ 

The  appropriate  value  of  Cf  for  use  in  this  solutioi^  is  obtained  from  thi 
relation  (see  Section  6): 

2l  =  1  C28) 

C-  1  +  0.633  C 

lO 

* 

The  dependence  of  6  and  on  mass  addition  rate  is  deduced  from  the 

approximation: 


(29) 


The  VIZAAD  results  are  compared  with  test  data  in  Figures  43  through  45 
(solid  lines).  The  method  introduces  a  dependence  on  bluntness  ratio. 
Agreement  with  measured  data  appears  to  be  better  for  the  sharp  cone  than 
for  Rif/Rg  =  0.077. 

The  dependence  on  mass  addition  rate  is  introduced  by  McDonnell-Douglas^22) 
through  the  following  relations: 


a* 

6* 

o 


+  c 


(30) 
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where 


~^—  -  1  +  0,68  ^  +  0.08  W®‘  ^ 


(31) 


and 


w  =  y.  .  /y^  , 

mj.  b.  1. 


In  the  present  case  the  value  of  the  ratio  W  was  unity.  Agreement  with 
test  data  is  good  for  the  mass  addition  rates  shown. 
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SECTION  VI 


SKIN  FRICTION  REDUCTION  CORRELATION  -  LAMINAR  FLOW 


Sufficient  data  to  assess  the  effect  of  mass  injection  on  the  skin  friction 
reduction  were  not  obtained  in  the  experimental  investigations(2b, 28)  discussed 
in  the  body  of  this  report.  However,  the  laminar  force  data  of  Ref.  24  in  con¬ 
junction  with  the  large  body  of  data  of  King(37)  provide  a  suitable  basis  for 
reassessing  the  correlations  which  currently  exist.  Since  the  skin  friction 
was  not  measured  directly  in  the  experiments,  it  was  obtained  indirectly  by 
calculation  as  follows 


'D 


TOTAL 


+  C 


(32) 


All  drag  coefficients  are  referenced  to  the  cone  base  area.  The  total 
drag  is  the  experimental  value  registered  by  the  drag  balance. 


The  pressure  drag  is  defined  as  the  streamwise  component  of  force  due  to 
the  slant-face  pressure  in  excess  of  free  stream  pressure.  This  can  be  ex¬ 
pressed  as: 


2  IT  tan^  6^  L 

1/2 \  J 


(P-PJ  xdx 


(33) 


The  pressure,  p,  referred  to  here  is  the  experimentally  determined  value. 


The  base  drag  corresponds  to  that  vdiich  occurred  in  the  experiment  as  a 
result  of  a  particular  method  of  mounting  the  models  in  the  wind  tunnel,  and 
consequently  has  no  significance  beyond  that  of  an  experimental  correction. 
Because  the  pressure  existing  on  the  base  and  inside  the  hollow  afterbody 
was  different  from  P®  ,  a  net  force  occurred  on  the  model.  The  base  drag 
coefficient  can  then  be  expressed  as 


1/2  P  A 


(34) 


The  negative  sign  is  inserted  because  pressure  on  the  base  actually  produces 
a  thrust  on  the  model. 


The  passage  of  injected  gas  from  the  model  surface  into  the  boundary  layer 
produces  a  drag  force  increment  which  is  small,  but  not  negligible.  If  we 
assume  the  validity  of  the  no-slip  boundary  condition,  then  the  injection 
velocity  must  be  normal  to  the  cone  surface.  The  injection  drag  coefficient 
component  may  be  computed  from  the  flux  of  momentum  across  the  porous  surface 
as 
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Shown  in  Figure  46  are  the  data  of  King(37)  plotted  versus  the  blowing 
parameter,  C  .  Data  were  obtained  for  two  values  of  the  free  stream  Mach 
number  (i.e.  =  3.93  and  5.64  or  Mg  “  3.75,  5.29)  and  for  two  injectant 

gases  (i.e.  Helium  and  air).  The  data  for  air  injection  show  no  apparent 
effect  of  free  stream  Mach  .nixnber  as  is  evident  in  the  Helium  data.  In 
addition,  it  appears  that  Helium  is  not  as  effective  as  air  in  reducing  the 
drag,  which  is  in  agreement  with  the  conclusions  of  Little  and  Griffith . 
Also  shown  in  Figure  46  is  the  data  sinnmary  of  Leadon  et  al^^^)  Kent,  et  al(28) 
and  some  unpublished  GE  data.  The  experimental  results  have  been  empirically 
curve-fitted  using  the  equation  for  Cf/Cf^  derived  in  Dorrance(^),  where 


expanding  and  retaining  only  linear  terms  in  C  »  there  results: 

1  (39) 

So  ' 


where  K  is  the  empirical  constant  that  achieves  a  best  fit  with  the  data. 
It  was  concluded  chat 

^  1  (^) 
CjJ*  "  1  +0.663  C 

is  a  "best-fit”  as  shown  in  Figure  46. 

The  empirical  techniques  of  Timner  et  al^^^  and  Costello^^^^  are  also 
shown  in  Fig.  46.  It  is  evident  that  both  of  these  methods  tend  to  under¬ 
predict  Che  skin  friction  drag  of  an  ablating  vehicle. 
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SUMMARY  AND  CONCLUSIONS 


Measurements  of  pitot  pressure  and  stagnation  temperature  profiles  as 
well  as  surface  pressures  and  temperatures  obtained  at  AEDC  under  1969 
STREET-G  and  1968  Mark  12  p/A  contracts  have  been  reduced  and  correlated. 
Values  of  boundary  layer  thickness,  displacement  and  momentum  thicknesses, 
form  factor,  profile  exponents  n  and  c  ,  and  location  of  boundary  layer 
transition  have  been  deduced  for  turbulent  boundary  layers.  In  addition, 
the  first  three  of  these  parameters  have  been  deduced  in  several  cases 
where  boundary  layer  flow  remained  laminar. 

The  data  treated  in  this  study  was  obtained  on  a  7.25  degree  sphere 
cone  with  bluntness  ratios  of  zero,  0.077  and  0.180.  Free  stream  Mach 
number  was  8.0  and  Reynolds  number  ranged  from  0.65  x  10^  to  3.7  x  10® 
per  foot.  Integrated  mass  addition  ratio  varied  from  zero  to  0.05. 

Results  deduced  from  experiment  are  compared  with  empirical  prediction 
methods  used  at  General  Electric,  McDonnell-Douglas  and  Philco  Ford.  All 
three  methods  appear  to  overpredict  the  effect  of  mass  addition  on  boundary 
layer  thickness  in  turbulent  flow.  A*new  set  of  equations,  accounting  for 
the  effect  of  mass  addition  on  6 ,  6  ,  and  6  is  proposed  for  turbulent 
boundary  layers.  In  addition,  a  correction  is  proposed  to  a  constant  used 
with  the  GE  thickness  predictions  at  zero  mass  addition.  The  proposed 
changes  result  in  an  improved  drag  and  plasma  attenuation  prediction  for 
flight  vehicles. 

In  laminar  bomidary  layers, currently  used  correlations  appear  adequate 
at  low  mass  addition  rates.  The  quantity  of  available  measured  data  was  not 
adequate  to  permit  an  evaluation  of  currently  used  correlations  at  high  mass 
addition  rates. 

In  several  turbulent  cases  measurements  were  compared  with  results  com¬ 
puted  with  the  GE  Equilibrium  Non-Similar  Boundary  Layer  Program  (ENSBL). 

a  finite  difference  solution  of  the  boundary  layer  equations.  Three  eddy 
viscosity  models  were  used  with  ENSBL,  a^  results  computed  with  each  are 
presented.  Good  agreement  is  noted  in  6  and  6*,  but  6  is  slightly  under¬ 
predicted  by  ENSBL. 
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TABLE  1.  KEY  FOR  BOUNDARY  LAiER  PROFILE  DATA 
A.  TURBULENT 
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Figure  2.  Schemalic  Drawing  of  the  Model  Configuration 
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Figure  4.  Displacement  Thickness  and  Mass  Addition  Rate  (ENSBL  Calculation) 
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Figure  9.  Transition  Detection  from  Pitot  Pressure  Survey  on  Model  Surface 
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Figure  11.  Effect  of  Mass  Addition  on  the  Boundary  Layer  Thickness  of  a 
Sharp  7. 25°  Half  Angle  Cone  (Re  =  3. 7  x  10%-!) 
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Figure  13. 
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Figure  14.  Effect  of  Mass  Addition  on  the  Boundary  Layer  Thickness  of  a 
Sharp  7.  25°  Half  Angle  Cone  (Re  =  2. 7  x  10®  ft~^) 
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Figure  18.  Effect  of  Mass  Addition  on  Velocity  Piofile  Exponent 
(Re«  =  2.7  X  10®  ft“^) 
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Figure  19.  Variation  of  Velocity  Profile  Exponent  with  Normalized 
Length  from  End  of  Transition 
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Figure  27.  Comparison  of  Total  Enthalpy  Profile  with  Crocco  Theory 
(Effect  of  Mass  Injection) 
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Figure  28.  Effect  of  Mass  Injection  on  Total  Enthalpy  Profile 
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Figure  31.  Comparison  of  Theory  with  Experimental  Viscous  Layer 

Thicknesses  (Case  B) 
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Figure  32.  Comparison  of  Theory  with  Experimental  Viscous  Layer 

Ihicknesses  (Case  C) 
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Figure  33.  Comparison  of  TTieory  svith  Experimental  Viscous  Layer 

Thicknesses  (Case  D) 
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Figure  35.  Slender  (R/V)  S  Band  Plasma  Attenuation 
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Figure  .“16.  Comparison  of  ENSBL  Program  Results  with  Experimental  Data 
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Figure  38.  Comparison  of  Predicted  Enthalpy  Profiles  with  Crocco 
Theory  and  Data  -  Zero  Mass  Injection 
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Figure  42.  Comparison  of  Theory  wift  Experimental  Viscous  Layer 

Thicknesses  (Laminar) 
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Figure  46.  Effect  of  Mass  Injection  on  Skin  Friction  Reductioii  -  Laminar 


AFPENDLX  A 

TABULATED  SUMMARY  OF  PROF  H  E  DATA 
BOUNDARY  LAYER  DATA  DEDUCED  FRW  PROFILES 


®c  “  7.25® 

N2  “  Injectant 


(P  “)e 

_ ^  1-T^SEC 

_ ^  FT^SEC 


Cfo(THKOR.)  X  103 
if.  2  X  (THEOR.) 

pCrp _ 

;  60  (INCHES) 


8*  (INCHES) 
00  (INCHES) 


j  H„=»oy®o 


;  n  (vflocity  profile 
O  exponent) 


«/«0 

■  *7  * 

I  » / « ■ 

!  «/«o 


CASE  A 

AXIAL  STATION  (INCHES) 

14.61  21,61  30.61  36.61 

6.63  9.24  13,3  16.0  18.2 

18.3  18.3  18.3  18.3 

_0 _ 0 _ 0 _ ^ 

_P _ 0 _ 0 _ 0_ 

.9561  .8949  .8415  .8133 


8.9  10.8 


.173  .270  .368  .415 

.0554  .0832  .148  ^8  .194 

.00529  .00856  .0119  _ jOl 

10.47  9.72  12.44  12.80 


7.2  6.65 


m 

u 


!  11/ Ho 

L- - 

!  .n 

j  ‘w/reo 

!  ‘Vi’* 


7.85 

7.62 

7.48 

1 

2.80 

3.00 

2.80 

1 

I 

!  2.8 

9 


BOUNDARY  LAYER  DATA  DEDUCED  FR^M  PROFILES 

CASE  B 


©c  «  7.25® 

N2  “  Injeccant 


Res 

X  10 

(/0u)_ 

LBM 

''0  ft^sec 

(Pv) 

LBM 

_ ^  FT^SEC 

^  =(pv)^ 

Cfo(THKOR.)  X  103 

C=— L  A  (THSOR.) 

Cfo _ 

6  (INCHES) 

6  *  (INCHES) 

0  (INCHES) 

H=  aye 

n  (velocity  profile 
exponent) 


_ AXUL  STATION  (INCHES)  _ 

14.61  21.61  30.61  36.61  AT  .^.1 

6.63  9.24  13.3  16.0  18,2 

18.3  18.3  18.3  18.3  18.3 

.00612  .00612  .00612  .00612  .00612 

3.344  3.344  3.344  3.344  3.344 

.9561  .8949  .8415  .8133  .7935 

.700  .747  .795  .822  .843 

.•216  .340  .440  .510  .515 

.0743  .112  .182  .223  .238 

.00701  .0112  .0156  .0183  .0191 

10.60  10.00  11.67  12.19  12.46 

9.25  7.65  6.15  5.75  6.60 


«/«o 

1.25 

1,26 

1.195 

1.23 

1.132 

5*/  6*0 

1.341 

1.346 

1.230 

1.253 

1.227 

1.325 

1.308 

1.311 

1.317 

1.282 

h/ho 

1.014 

1.03 

.938 

.952 

.957 

n/no 

1.04 

.708 

.854 

.865 

.924 

Ty/TOf, 

7.44 

7.12 

6.85 

6.75 

6.73 

2.91 


2.92 


2.90 


2.90 


2.80 


Injectant 


..  nn-6 


(O  u)^  LBM 
6o  .  ' 


Cfo(THEOR.)  X  103 


.  2  X  (THEOR.) 

'  Cfo 


(INCHES) 


6*  (INCHES) 


0  (INCHES) 


H=  6*/6 


n  (velocity  profile 
exponent) 


*/  * 

6  /  6  o 


8/  9o 


h/Ho 


n/no 


VTeo 


P„/P, 


14.61 

21.61 

30.61 

36.61 

6.63 

9.24 

m 

16.0 

18.3 

18.3 

18.3 

18.3 

,01115 

1  .01115 

.01115 

.01115 

6.093 

5.093 

6.093 

6.093 

,9561 

.8949 

1 

.8415 

.8133 

-.274 

1.362 

1.448 

1.498 

.230 

.378 

.505 

.558 

1085 

.137 

.218 

.251 

,00866 

.0142 

.0193 

.0222 

12.53 

9.65 

11.30 

11.30 

11.9 

6.65 

5.35 

4.95 

1.33 

1.40 

1.372 

1.345 

j 

1.958 

1.647 

1.473 

1.410 

1.637 

1.659 

i;622 

1.600 

1.20 

.993 

.908 

.883 

.39 

.615 

.743 

.745 

.02 

6.70 

6.42 

6.31 

8 


M  -  8 

m 

8c  -  7.25® 

No  -  Injectant 


BOUNDARY  LAYER  DATA  DEDUCED  FROM  PROFILES 


CASE  D 

ant  AXIAL  STATION  (INCHES) 


14.61 

21.61 

30.61 

36.61 

41.51 

Reg  X  10‘^ 

6.63 

9.24 

13.3 

16.0 

18,2 

(Ou)  LBM 

o  FT^SEC 

18.3 

18.3 

18.3 

18^3 

18.3 

(Pv) 

^  FT^SEC _ 

.0303 

.0303 

.0303 

.0303 

1 

.0303 

X  .(Pv)^/«>u)^^  X  10^ 

16.56 

16.56 

16.56 

16.56 

16.56 

Cfo(THF,OR.)  X  103 

.9561 

.8949 

.8415 

.8133 

.7935 

r_  2  X  rTHEOR.") 

‘  Cfo 

3.463 

3.700 

3.935 

4.072 

4.1/3 

1 

6  (INCHES) 

.650 

.775 

.885 

6  *  (INCHES) 

.209 

.365 

.425 

f)  (INCHES) 
H-  6*1$ 

n  (velocity  profile 
exponent) 


H/io 


*/  * 

«/  6o 


®/«o 


h/ho 


n/ Hq 


T^/Tcq 


,0219  .0355  .0393 
9.54  10.28  10.8 
4.05  3.6  4.1 


1.766  1.87  1.945 


1.412  2.051  2.190 


1.840  2.554  2.640 


3 


.563 

.541 

1 

.573 

5.12 

5.01 

4.56 

2.86 

2.82 

2.74 

BOUNDARY  LAYER  DATA  DEDTCBD  EROH  PROFILES 


N2  “  Injectant 


(P  u)  LBM 
60“"^ 


X  MP%  //)u)e^  X  10^ 


Cfo(THKOR.)  X  103 


(THEOR.) 


(INCHES) 


(INCHES) 


0  (INCHES) 


1=  s*/e 


n  (velocity  profile 
exponent) 


«/*» 


«7 


®/«0 


h/ho 


T„/Tej 


CASE  E 

AXIAL  STATION  (I^HES) 
14.61  21.61  30.61  36.61 


6.63  9.24  13.3  16.0  18.2 


18.3  18.3  18.3  18.3  18.3 


.060  .060  .060  I  .060  |  .060 


32.79  32.79  32,79  32.79  32.79 


6.858  7.327  I  7.792  |  8.063  I  8.264 


.925  1.175  1.32 


.535  .649 


.0492  .0599  .0660 


10.87 


11.17 


2.70  2.50  3.10 


2.51  2.835  2.90 


3.615  3.646  3.799 


4.134  4.309 


.874  .846 


.375  .376  .434 


4.04  I  3.99  3.94 


3.10  3.10  2.71 


M*  -  8 

Oq  -  7.25° 

N2  “  Injectant 

"TTT” 


(0  u)  LBM 
®0  FT^SEC 


(P  V)  LBM 
_ ^  FT^SEC 

A  =  (pu  /(Pu).  X  1( 
_ ’  *^o 

Cfo(THK0R.)  X  10^ 

^_2 _ ^A  (THEOR.) 

Cfo _ 

80  (INCHES) 

5*  (INCHES) 

9q  (INCHES) 

H  =6*/9 
0  0/  o 

n  (velocity  profile 
o  exponent) 


*/  * 

6/  «o 


^ho 


h/ho 


n/oo 


Ty/TCp 


Pw/P« 


BOUNDARY  LAYER  DATA  DEDUCED  FROM  PROFILES 


CASE  F 

ant  _  AXIAL  STATION  (INCHES) 


21.61 

30.61 

36.61 

41.51 

7.23 

10.25 

12.25 

13.89 

13.21 

13.21 

:  13.21 

1 

i 

1  13.21 

0 

0 

0 

0 

0 

0 

0 

0 

.9374 

.SS42 

.8568 

.8375 

0 

0 

0 

0 

.208 

.357 

.393 

.458 

.0698 

.137 

.158 

.195 

.00593 

.0112 

.0131 

.0153 

! 

11.77 

12.23 

12.06 

12.74 

6.7  7.9  7.2  I  7.0 


7.93 

7.72  ! 

7.65 

7.61 

2.85 

2.86 

2.85 

2.86 

BOUNDARY  LAYER  DATA  DEDUCED  FROM  PROFILES 


BOUNDAKf  LAYER  DATA  DEDUCED  FROM  PROFILES 


M  -  8  CASE  G 


8c  -  7.25® 


WC  • 

K2  -  Inject*nt 

AXIAL  STATION  (INCHES) 

14.61 

21.61 

30.61 

36.61 

41.51 

Reg  X  10”^ 

7.23 

10.25 

12.25 

13.89 

(Ou)  LBM_ 

o  FT^SSC 

13.21 

13.21 

13.21 

13.21 

(.rJv)  L8M_ 

^  FT^SEC 

.00379 

.00379 

.00379 

.00379 

2.869 

2.869 

-  2.869 

2.869 

Cfo(THEOR.)  X  103 

.8842 

.8375 

5,  2  (THEOR.) 

Cfo 

• 

.649 

.685 

fi  (INCHED) 

.400 

r  '  —  ■ 

.505 

d  *  (IKCHES) 

.169 

. 

.224 

®  (INCHES) 

«» 

.0140 

.0186 

H=  b*/B 

12.05 

12.04 

n  (velocity  profile 
exponent) 

• 

6.9 

7.3 

»/*o 

1.12 

1.105 

6y  6*0 

1.233 

1.149 

«/9.o 

1.250 

1.216 

h/ho 

.985 

.945 

n/no 

.873 

1.04 

7.33 

7.24 

Vp. 

2.83 

2.90 

BOUNDARY  LAYER  DATA  DEDtXJED  FROM  PROFILES 


M  -  8  CASE  H 

m  ““ 

Be  -  7.25®  _ 


^2  •  Injeettnt 

AXIAL  STATION  (INCHES) 

14.61 

21.61 

30.61 

36.61 

41.51 

Re.  X  10“^ 

7.23 

10.25 

12.25 

13.89 

(PU)  LBM 

_  O  FT^^SEC  _ 

13.21 

13.21 

13.21 

13.21 

.00803 

.00803 

.00803 

.00803 

Xs(pv)^  ^0^ 

6.079 

6.079 

6.079 

6.079 

Cfo(THEOR.)  X  103 

.9374 

.8842 

.8568 

.8375 

X  (raeoR.) 

Cfo 

1.297 

1.375 

1.419 

1.452 

6  (INCHES) 

.338 

.500 

.570 

.658 

6  *  (INCHES) 

.128 

.213 

.245 

.282 

®  (INCHES) 

.0118 

.0175 

.0211 

.0243 

a=  eye 

10.85 

12.17 

11.61 

11.60 

n  (velocity  profile 
exponent) 

n 

n 

6.2 

6.45 

»/'o 

1.425 

1.40 

1.45 

1.445 

«*/  6*0 

1.834 

1.555 

1.551 

1.446 

»/«•. 

1.990 

1.563 

1.611 

1.588 

h/ho 

.922 

.945 

.963 

.910 

n/nj, 

.955 

,898 

IliKHIIl 

m 

VTeo 

7.07 

6,91 

6.84 

6.79 

VP, 

2.90 

2.80 

2.84 

2.68 

88 


BOUNPARY  lAYER  DATA  DEDUCED  FROM  PROFILES 


VTeo 


5.71 


5.66 


I 


BOUNDARY  LAYER  DATA  DEDUCED  FROM  PROFILES 


M  "  8  CASE  M 


Sc  -  7,25® 


N2  *  Injectnnt 

AXUL  STATION  (INCHES) 

14.61 

21.61 

30.61 

36.61 

41. SI 

Re,  «  10‘‘ 

7.23 

10.25 

12.25 

13.89 

(Pay  LBM 

®oft^sec 

13.21 

13.21 

13.21 

13.21 

(Ov)  LBM 

_ .y.JT^SEC _ _ 

.00422 

.00422 

.01265 

.01265 

3.195 

3.J95 

9.576 

9.576 

Cfo(THEOR.)  X  103 

.9374 

.8842 

.8568 

.8375 

X  (THEOR.) 

Cfo 

.682 

.723 

2.235 

2.287 

6  (INCHES) 

.290 

.445 

.505 

.590 

6*  (INCHES) 

.0943 

.175 

.232 

.281 

9  (INCHES) 

.00978 

.0153 

.0199 

.0237 

H=  5*/© 

9.64 

11.44 

11.66 

11.86 

n  (velocity  profile 
exponent) 

6.8 

7.15 

7.10 

4.95 

5.30 

»/«o 

1.395 

1.245 

1.285 

1.29 

6*/  »*c 

1.350 

1.277 

1.468 

1.441 

1.650 

1,366 

1.519 

1.549 

h/Ho 

.819 

.935 

■i 

■1 

n/no 

QQQH 

.900 

,.688 

.757 

Ty/Tco 

7.47 

7.30 

6.37 

6.33 

BOUNDARY  UYER  DATA  DEDXED  FROM  PROFILES 


M  -  3  CASE  N 

•  — 

Cc  -  7,25°  _ _ 


V 

N2  ~  Injectant 

AXIAL  STATION  (INCHES) 

14.61 

21.61 

30.61 

36.61 

wmm\ 

Reg  X  10“^ 

7.23 

10.25 

12.25 

13.89 

(pu)  LBM 

oft^sec. 

13.21 

13.21 

13,21 

13.21 

(Ov)  LBM 

w  ft-^sec 

.01158 

.01158 

.03474 

.03474 

Xs(pv)  /(0“)o  *  10^ 

IV  o 

8.766 

8.766 

26.30 

26.30 

Cfo(THEOR.)  X  10^ 

.9374 

.8842 

.8568 

,8375 

2  X  (THEOR.) 

Cfo 

1.983 

6.280 

6  (INCHES) 

m 

.608 

.925 

6  *  (INCHES) 

.240 

m 

.506 

®  (INCHES) 

• 

.0227 

m 

.0398 

a: 

It 

10.57 

12.71 

n  (velocity  profile 
exponent) 

■1 

■HI 

»/'o 

1.705 

2.02 

*/  * 

6/60 

1.752 

2.595 

9/9-0 

2.027 

2.601 

h/Ho 

.864 

.998 

n/no 

.810 

.450 

810 


450 


BOUNDARY  LAYER  DATA  DEDUCED  FROM  PROFILES 


M  -  8 

Oq  -  7.25® 

CO  -  injftctant 

A 


case  bb 

_ AXIAL  STATION  (INCHES) 

14.61  21.61  r  30.61  36.61 


BOUNDARY  LAYER  DATA  DEDUCED  FROM  PROFILES 


M  -  8  CASE  GG 

•  ■ 

0c  -  7.25°  _ 


N2  “  Injectant 

AXIAL  STATION  (INCHES) 

14.61 

21.61 

30.61 

36.61 

Riia 

Reg  X  10”^ 

11.98 

13.58 

(pu)  lbm 
_ o  FT^SEC _ 

12.92 

12.92 

DBIBHi 

.00614 

.00595 

X*(pv)  / 

•a  • 

4.752 

0 

CfjjCTHEOR.)  X  103 

.8568 

.8375 

2  ^  (THEOR.) 

Cfo 

1.110 

0 

6  (INCHES) 

.510 

.550 

t*  (INCHES) 

- 

.133 

®  (INCHES) 

.0120 

X 

II 

- 

11.07 

n  (velocity  profile 
exponent) 

- 

6.8 

»/'o 

1.200 

«y  6*0 

» ho 

h/ho 

n/n<, 

7.08 

7.03 

Pw/P. 

2.71 

2.61 

96 


aOUNDARY  LAYER  DATA  DEDUCED  FROM  PROFILES 


BOUNDARY 


$c  -  7.25® 

N2  “  lnject«nt 


Reg  X  10“^ 


(PuV  LBM 

■  ...rofXtSE^L, _ 

(pv)  ^bm 

. y-jFT.^saL. _ 

X2(pv)  /(OU)  X  10^ 

1  w  * 


Cfo(THEOR.)  X  10^ 


Cf. 


(THEOR.) 


d  (INCHES) 


fi  *  (INCHES) 


®  (INCHES) 


H=  aye 


n  (velocity  profile 
exponent) 


99 


Cfo(THEOR.)  X  103 


C»~J_x  (TOEOR.) 

Cfo 


(INCHES) 


*  (INCHES) 


®  (INCHES) 


H=  sye 


n  (velocity  profile 
exponent) 


«/'» 


*/  * 
6/60 


n/fiQ 


^w/Tco 


VP. 


7.86 

7.65 

2.22 

2.32 

BOUNDARY  UYER  DATA  DEDUCED  FROM  PROFILES 


-  8  CASE  SL 


M  -  8 

6c  -  7.25® 

N2  *  Injectant 


BOUNDARY  LAYER  DATA  DEDUCED  FROM  PROFILES 

CASE  TL 

> 

:8nt  AXIAL  STATION  (INCHES) 


(Pu)  LBM 
®o  ft2sec 


(pv)  LBM  . 

_ y-FT.^SEC _ 

Cfo(THEOR.)  X  103 

;^  .A  X  (THEOR.) 
Cfo 


6  X  '3) 
fi  awCHES) 

®  (INCHES) 

H=  6*/e 

n  (velocity  profile 
exponent) 


*/  * 

6  /  6  o 


«/«.„ 


h/h„ 


n/ny 


VTeo 


30.61 

36.61 

41.51 

2.619 

3.217 

3.55 

3.056 

3.056 

3.056 

.00116 

1 

1 

.00118 

0 

3.796 

3.861 

.. 

.711 

.650 

.610 

1.068 

1.188 

0 

.375 

.430 

.488 

.249 

.259 

.210 

.0141 

.0X80 

.0221 

17.63 

14.40 

9.46 

1.82 


.928 


2.62 


7.14 

7.65 

2.68 

2.58 

1 


00  -  7.25® 


N2  ^  InJectAnt 


BOUNDARY  LAYER  DATA  DEDUCED  FROM  PROFILES 


CASE  UL 

► 

:«nt  I  AXIAL  STATIW  (INCHES) 


14.61 

21.61 

30,61 

36. £ 

>1 

Re,  X  10“^ 

2.619 

3.21 

7 

3.55 

(PU).  LBM 

®0  PT^SEC 

3.056 

3.05 

6 

3.056 

(pv)  yw 

w  rT2sEC 

Xs(pv)  /<Pu* ,  X  10^ 

Cfo(THEOR.)  X  103 

C— i-  X  (tHEOR.) 

_ 

S  (INCHES) 

0  *  (INCHES) 

®  (INCHES) 

H=  aye 

n  (velocity  profile 
exponent) 


.00245  .0026/^  |  0 

8.017  8.639 

.711  .650  .610 

2.255  2,658  0 

.460  .50 

.371  ,342  I  .264 

.0141  .0198 

26.41  17.27  13.35 


1.99 

1.98 

1.95 

1.90 

1.61 

1.17 

1.94 

2,50 

2.J5 

.984 

.644 

.498 

5.20 

7.35 

7.14 

i 

2.91 

2.77 

2.53 

BOUNDARY  LAYER  DATA  DEDUCED  FROM  PROFILES 


M  -  8  CASE  VL 

m 

Be  -  7.25®  _ 


#  •  •.v 

N2  •  Injectant 

AXUL  STATION  (INCHES) 

14.61 

21.61 

30.61 

36.61 

41.51 

Reg  X  10"^ 

.842 

1.277 

(Ou)  LBM 

®0FT^SEC 

1.668 

1.952 

(flv)  LBM 

.  W  FT.2SEC 

0 

0 

0 

0 

Cfo(THEOR.)  X  103 

1.221 

.994 

2  ^  (THEOR.) 

Cfo 

0 

0 

6  (INCHES) 

.225 

.250 

B*  (INCHES) 

.150 

.172 

®  (INCHES) 

.0183 

.0163 

H=  6*fB 

8.20 

10.59 

1 


I 

\ 


n  (velocity  profile 
exponent) 


T^/Tcq 

VP. 


6.22 


6.22 


BOUNDARY  LAYER  DATA  DEDUCED  FROM  PROFILES 


CASE  XL 


Be  -  7.25® 


^2  “  InJectant 

AXIAL  STATION  (INCHES) 

14.61 

21.61 

30.61 

36.61 

Reg  X  10"^ 

- 

.842 

1.277 

(Pu)  LBM 

®  ^r'^SEC 

1.668 

1.952 

(Pv)  LBM 

w  ft^sec 

.00263 

.00210 

X  =  (pv)  /  (Ou)  X  10^ 

w  '  *=o 

15.77 

10.76 

Cfo(THEOR.)  X  103 

1.221 

.994 

C* _ Lx  (THEOR.) 

Cfo 

2.583 

2.165 

fi  (INCHES) 

.575 

.600 

B  *  (INCHES) 

.455 

.476 

®  (INCHES) 

.0342 

.0336 

H  =  6*^6 

13.31 

14.15 

n  (velocity  profile 
exponent) 

- 

«/'o 

2.56 

2.40 

+/  ★ 

6  /  6  o 

1  ' 

3.03 

2.77 

L  .  ■  _ 

»/«« 

1.87 

2.06 

h/Ho 

1.62 

1.34 

n/no 

- 

- 

Ty,/TCQ 

5.10 

4.90 

Pw/P. 

2.81 

2.48 

106 


BOUNDARY  LAYER  DATA  DEDUCED  FROM  PROFILES 


0c  -  7.25® 

N2  “  InJectAnt 


(pu)  lbm 

”0  FT^SEC 

(pv)  lbm 

_ ^  FT^SEC 

Aa(Pv)  /(0«)e  X  10^ 
w  ‘  _ 

Cfo(THEOR.)  X  103 

r«— L  X  (THEOR.) 

:_cfo_ _ 

0  (INCHES) 

0*  (INCHES) 

®  (INCHES) 

H=  6*/e 

n  (velocity  profile 
exponent) 


«/«o 


*/  * 

0  /  0  o 


»/ho 


"/"o 


VTeo 


VP. 


CASE  ~ 

m  g  0  Case  Corresponding  to  UL,  TL  (VIZAAD) 
_  AXIAL  STAflON  (INCHES) 


30.61 

36.61 

41.51 

2.66 

1  3.16 

3.60 

3.115 

1 

!  3.105 

3.105 

.701  .644  .604 

0 _ 0 _ 0 

.231  .252  .  269 

.195  .212  .226 

.00726  .00792  .00844 
26.84  26.83  26.82 


7.28  7.28  7.28 


